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rpl ' Abstract. We present formulas for computing the resultant of 

'^H , sparse polynomials as a quotient of two determinants, the denom- 

Osl ' inator being a minor of the numerator. These formulas extend the 

Cn , original formulation given by Macaulay for homogeneous polyno- 

mials. 
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1. Introduction 



Let ^0; • • • ; -^n be finite subsets of Z" and consider n + 1 polyno- 
mials /o, . . . , /n in n variables such that supp(/j) C Ai, i = 0, . . . , n. 
CN I The sparse resultant is an irreducible polynomial in the coefficients of 

/o, . . . , fn, which vanishes if the system fi = 0,i = 0, . . . ,n has a solu- 
00 ■ tion in an algebraically closed field. It will be denoted by Res^(/o, ... , /„ 

, where A := (A, • • • , An). 

Q I Resultants eliminate the input variables, so they are also called elim- 

inants. They have been used in the last decade as a computational tool 
for elimination of variables and for the study of complexity aspects of 
^-^ I polynomial system solving. This has renewed the interest in finding ex- 

^ ' plicit formulas for their computation (see [AS, Canl, Can2, CEl, CE2, 

CDS, CLO, DD, Emil, EM, KPS, Laz, Ren, Roj, Stul, Stu2, Stu3, 
ZCG]). 

'k>i I The study of resultants goes back to the classical work of Sylvester, 

^ ■ Bezout, Cay ley, Macaulay and Dixon in the context of homogeneous 

polynomials ( [Syl, Bez, Cay, Mac, Dix]). The sparse resultant, a gen- 
eralization of the classical one, first appeared in the study of hypergeo- 
metric functions and ^discriminants done by Gelfand, Kapranov and 
Zelevinski a few decades ago ([GKZl, GKZ2]). 

The first effective method for computing the sparse resultant was 
proposed by Sturmfels in [Stul]. In [CEl, CE2], Canny and Emiris 
gave algorithms for computing square Sylvester style matrices with de- 
terminants equal to non-zero multiples of the resultant. By a Sylvester 
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style matrix we mean the matrix in the monomial bases of a hnear map 
given by a formula as follows: 

{90,91, ■■■ ,9n) ^ Y.9ifi- 

Here, So, ■ ■ ■ Sn, S are finite sets of monomials in a ring of Laurent poly- 
nomials IK [xi,x^^, . . .x„,x~^] , and Sb denotes the K— vector space 
generated by B. This construction was generalized by Sturmfels in 
[Stu2] and it was pointed out in [Emil, CE2] that the extended formu- 
las, when applied to the classical case, give Macaulay's original formu- 
lation (see [Mac]). 

However, Macaulay succeeded in giving an explicit formula for the 
extraneous factor appearing in his own formulation, i.e. he showed 
that, in the classical case, the ratio 

det( Sylvester matrix) 
Res^(/o,--- ,fn) 

is a minor of the Sylvester matrix (see [Mac]). This was conjectured 
to happen in the sparse case, but no proof of it was available (see 
[CEl, CE2, CLO, Emil, EM, GKZ2, Stu2]). For instance, in [CLO, 
Chapter 7], we may read: 

One of the major unsolved problems concerning sparse resultants is 
whether they can he represented as a quotient of two determinants. In 
the multipolynomial case, this is true by Theorem (4-9) of Chapter 3. 
Does this have a sparse analog? Nobody knows! 

In [GKZ2, Introduction], the following is written: 
Macaulay made another intriguing contribution to the theory by giv- 
ing an ingenious refinement of the Cayley method [Mac]. It would be 
interesting to put his approach in the general framework of this book. 

In [Stu2, Corollary 3.1], we also find: 
It is an important open problem to find a more explicit formula for 
Plo,5^ m the general case. Does there exist such a formula in terms of 
some smaller resultants? 

The main contribution of this paper is a positive answer to this 
question, i.e. a generalization of Macaulay's classical formulas to the 
sparse case by means of an explicit algorithm which produces square 
Sylvester style matrices. The determinant of each of these matrices 
is a non trivial multiple of Res^(/o, . . . ,/n)- Moreover, we succeed 
in describing the extraneous factor of our formulation (i.e. the ratio 
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which appears in (1)) which again happens to be the determinant of a 
submatrix of the Sylvester matrix. 

The paper is organized as follows: in Section 2, some notation and 
preliminaries are introduced. In Section 3, we explicitly construct 
Sylvester style matrices for generalized unmixed families of polyno- 
mials and prove that our algorithm produces formulas "a la Macaulay" 
for computing the sparse resultant in this case. 

Section 3.2 deals with the general case, and may be regarded as an 
extension of the previous section. The algorithms are illustrated with 
examples at the end of both sections. 

Acknowledgements: I am grateful to Alicia Dickenstein who brought 
my attention to this problem and to loannis Emiris for helpful com- 
ments. I also wish to express my deep gratitude to David Cox for 
for his thorough reading of preliminary drafts of this paper and very 
thoughtful suggestions for improvement. 

This research began during the Long Semester Program in Symbolic 
Computation in Geometry and Analysis held at MSRI in the Fall Se- 
mester of 1998. I am grateful to the organizers for their help and 
support. I am especially grateful to Bernd Sturmfels for helpful con- 
versations during those days. 

2. Preliminaries 

We review here some definitions and properties of convex polytopes 
and sparse resultants. More details and proofs can be found in [CLO, 
EM, GKZ2, Stul, Stu2]. 

Let ^0) • • • ) -^n be finite subsets of the lattice Z". Set rrii := ^{Ai), 
m := X]r=o"^* ^"^^ Qi '■~ conv{Ai), i = 0, . . . ,n. Here, conv{') denote 
convex hull in C]^ := £ M, where C is the afiine lattice generated by 
V A- 

For any subset J C {0, 1, . . . , n}, consider the affine lattice generated 
by ^j^jAj, and let rk{J^ be the rank of this lattice. For every a G Ai^ 
we shall introduce a parameter Ci^a- Consider the family of generic 
polynomials: 

(2) fi (xi, . . . , x^) = ^ Ci^a x" (i = 0, . . . , n) . 

Let IK be an algebraically closed field. The vector of coefficients {ci^aJa^A 
of such a family defines a point in the product of K— projective spaces 
P;^""^ X ... X Pk""^ Let Z denote the subset of those families (2) which 
have a solution in (K*)" . Here, IK* denotes the torus K \ {0}. Finally, 
denote by Z the Zariski closure of Z in P^""^ x . . . x Pk""^ 
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Theorem 2.1. [GKZ2, Stu2] The projective variety Z is irreducible 
and defined over Q. Its codimension in P^"^ x . . . x P^" "^ equals the 
maximum of^{I)—rk{I), where I runs over all subsets o/{0, 1, . . . ,n}. 
The variety Z has codimension 1 if and only if there exists a unique 
family {Ai}i(=i such that 

1. rk{I) = #(/) - 1, 

2. rk{J) > i^{J), for each proper subset J of I. 

Definition 2.2. [Stu2] If / satisfies the conditions 1 and 2 given in the 
previous theorem, then the family {Ajjjg/ is said to be essential. 

Note that if each Qi is n— dimensional, it is easy to check that the 
unique set satisfying both conditions is / = {0, 1, . . . ,n}, so in this 
case Z has codimension 1. 

The sparse mixed resultant Res^(/o, . . . , fn) is defined as follows: 
if codim{Z) = 1, then Res^(/o, . . . , /„) is the unique (up to sign) 
irreducible polynomial in l^lci^a] which vanishes on Z. If codim{Z) > 2, 
then Res_4(/o, . . . , /„) is defined to be the constant 1. 

Theorem 2.3. [PS, Stu2] // the family of supports {Aq, . . . ,An} is 
essential, then for i = 0, . . . ,n, the degree of Res_4(/o, . . . , /„) in the 
coefficients of fi is equal to the normalized mixed volume 

MV{Qo, ... , Qi-i, Qi+i, ... ,Qn) := 



vol{V) 
where vol{) stands for the euclidean volume in the real vector space 
Cr, and V is a fundamental lattice parallelotope in C. In general, if 
there exists a (unique) subset {Ai}i(zi which is essential, the sparse 
mixed resultant coincides with the resultant of the family {fi : i G /}, 
considered with respect to the lattice Eie/'^*- 

Example 2.4. Set Aq = Ai = A2 = {(0, 0), (0, 1), (1, 0), (1, 1)}, and 
consider the family 

/o = Co,(o,o) + Co,(i,o)2^i + Co, (0,1)3^2 + Co,(i,i)a:i X2, 
/i = Ci,(o,o) + Ci,(i,o)a;i + Ci,(o,i)X2 + Ci,(i,i)Xi X2, 

/2 = C2,(0,0) + C2,(l,0)Xi + C2,{0,1)X2 + C2,(l,l)Xi Xg. 

Here, the Newton polytopes Qq, Qi, Q2 are equal to the unit square 
C := [0, 1] X [0, 1] whose vertices are precisely the points in the common 
support; and we have that 

degcoeff /, (Res^(/o, A, A)) = MV{C, C) = 2,t = 0, 1, 2. 
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A nice formula due to Dixon ([Dix]) allows us the computation of the 
resultant as follows: 

/ Co,(o,0) Co.d.o) Co.fo.l) co.a.ii \ 

Cl,{0,0) 



Res^(/o,/i,/2) = det 



v 



^2,(0,0) 







co,(i,o) 

Cl,(l,0) 
•^2,(1,0) 
Co, (0,0) 
Cl,(0,0) 
C2,(0,0) 



co,(o,i) 

'^2,(0,1) 







•^0,(1,1) 

^1,(1,1) 

^2,(1,1) 
C0,{0,1) 

"^2,(0,1) 







C0,(1,0) 
•^1,(1,0) 
"^2,(1,0) 







co,(i,i) 

^2,(1,1) 



/ 



Example 2.5. Let 



A = {(0,0), (2,2), (1,3)}, 
A = {(0,0), (2,0), (1,2)}, . 
A2 = {(3,0), (1,1)}. 



Consider the family 



/o 
/i 

/2 



«! + a2xl X2 + Q;3a;i X2, 



Pi + p2xi + [3^xi x\, 

1\X\ + l2X\ X2- 

A straightforward computation shows that 



MV[Q^, go = 7, MV{Q^, Q2) = 7, MV{Q^, Q2) = 5, 

and the sparse resultant equals 

«?/5l7h2 + 3ata2m^hl + SalalP^^t^l 
-13ala2asP!f32f3hhl ' M(4PiPWhhl + 6a?ai/3?/32/9|7i7l 
+alalPlP,^l - alalasPfPhhl + 5a?a2«i/5i^/5|7h2 
-ala2alP,PlP,jhl + UajalPfPlP^l^l + ajalPhnl 
-2a,aiPlPlMhl - ^a,ala,PlPhhl + al^f^.^hl 
+2aialal(3f(3lf3s^f^l - 2aia2al(3f(3hhl - 7aiaif3f(32(3s^h2 
+alalPlP2lhl + ^Whl 

For an explicit computation of this resultant, see [Stu2]. 

As usual, to define a face of a polytope Q C £r C M", let w be a 
vector in M". Set 



mqiv) 



min{(g,t;)}. 



and call 

Qv -■= Q ^ {m e M" : (m, v) = mQ{v)} 

the face of Q determined by v. The vector v will be called an inward 
normal vector of Qy. If dim.{Qy) = n — 1, then Q^ will be called facet. 
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3. Macaulay style formulas for generalized unmixed 
families of polynomials 

In this section, we will construct square Sylvester style matrices in 
the case all Qi are integer multiples of a fixed polytope P, i.e. there 
exists positive integer numbers /cq, A;i, . . . , A;„ such that 

Qi = conv{Ai) = kiP, i = 0, . . . ,n. 

This case is treated in [CDS] in a toric setting. These families may 
be identified with homogeneous polynomials in the coordinate ring 
of a projective toric variety (in the sense of [Cox]), with "degrees" 
ttO) • • • , «n, where every a, is Q— ample (see [Ful]). We shall call them 
generalized unmixed polynomials because they contain the well-known 
unmixed family of polynomials, which is the case when all input sup- 
ports Ai coincide. 

The matrices to be constructed here also generalize the formulas 
given by Macaulay in [Mac] in the homogeneous case, where all the 
supports are multiples of the standard simplex 

n 

Sn := {{qi, . . . ,g„)GM": < g, < 1, ^ g. < 1}. 

Warning: Some care must be taken with taking convex hulls, because 
the sparse resultant depends strongly on the finite data {Aq, Ai, . . . , An) 
and different families of input supports may give the same polytopes 
Qo,--- ,Qn (see [Stu2]). 

Remark 3.1. It is easy to see that, in the generalized unmixed case, 
because of Theorem 2.1, in order to have a nontrivial resultant, P 
must be n— dimensional. It is also clear that, in this case, C^ = IR"- 

Given A G Q>o and a generic 5 e C (g) Q as in [CEl, CE2, Stu2], 
our algorithm will produce a Sylvester style matrix M whose rows and 
columns will be indexed by the integer points in 

(3) S:={{ko + k, + ... + kn + X)P + 6)nC, 

and whose determinant will be a nonzero multiple of the sparse re- 
sultant. Moreover, we shall be able to identify the extraneous factor 
det(M)/Res^(/o, . . . , /„) as a minor of this determinant. 

The algorithm is recursive in the dimension of the polytope P, and in 
its intermediate steps, uses the mixed-subdivision algorithm of Canny 
and Emiris (see, [CEl, CE2, Stu2]) in order to refine the subdivision 
(see the comments in Remark 3.15). 



^0 -^ 

bo ^ 

b ^ 


R 
1 



if 6 ^ bo 


b ^ 


M 



ybeAi, i = l,...,n 


ViQ) - 
b ^ 


M 



V6 G V{Q). 


Ji{b)bGA,,i= 


0,1,... 


,n, C.(6),,v(Q) ) G M™+#^(«) 
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Remark 3.2. It was stated in the introduction that S should be a set 
of monomials. Indeed this is true provided that we identify an integer 
point a & S with the Laurent monomial x°'. 

3.1. Constructing the Matrix M. Given A and S as before, set 
Q := XP. Let V{Q) C £ ® Q be the set of vertices of Q. Observe that 
they are not necessarily integer points. 

Choose a vertex bo G ^o? and consider the following lifting functions: 

UJo 



(4) u. 



u 



Set 

and consider the lifted polytopes in M"+^ : 

Qi,n ■= conv{{a,uJi{a)) : a e Ai} 
Qn := conv{{b,u{b)) : beV{Q)}. 

By projecting the upper envelope of Qj,^ (resp. Qn) we get a coherent 
mixed decomposition Aj^n (resp. Af^) of the polytopes Qi (resp. Q). 
The cells in this decomposition are the projections of precisely those 
faces of Qi^n and Qn on which a linear functional with negative last 
coordinate is minimized (see [Stu2]). 

Similarly, by projecting the upper envelope of the Minkowski sum 

Qo,n + Qi,n + • • • + Qn,n + Qn, 
we get a coherent mixed decomposition of 
(5) Q:=Qo + Qi + ... + Qn + Q. 

Each cell in this decomposition is of the form 

F = Fo + Fi + . . . + F„ + F, 

where Fj (resp. F) is a cell in Aj^^^, (resp. A^^). 

Because the lifting functions we have used in (4) are most of them 
trivial, we can characterize all the proper cells in Ai q, . . . , A„ q, Aq : 

• the whole polytopes Qi, ■ ■ ■ ,Qn,Q corresponding to the linear 
functional associated to the vector (0, —1) G M" x M; 
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• for every v E W^ \ {0}, the faces Qi^, • • • , Qnv^ Qv associated to 
the vector (f,Q;), where a is any negative number. 

On the other hand, on Ao,n, the cell corresponding to (0, —1) is the 
singleton {bo}, and it is easy to check that every cell of dimension n 
in its decomposition is the convex hull of bo and a facet of Qq not 
containing this point. We shall call it Fo,t,, where v denotes the integer 
primitive inward normal of the facet. 

So, we can characterize all maximal (i.e. of dimension n) cells in the 
polyhedral decomposition of Q as follows: 

• {bo} + Qi + . . . + Qn + Q- This shall be called the primary cell. 

• Fo,v + Qiv +• • •+ Qnv + Qv, for some v G M". These will be called 
secondary cells, and are associated to a nonzero vector f G M". 

The following Lemma will be useful in the sequel: 

Lemma 3.3. Suppose that dim(Fo,t,) = n. Then, for every q E Q, if 
Fq^v n {m G Cm. '■ {iTT', v) = q\ ^ 0, there is \q G Q>o such that the 
intersection is a polytope congruent to XqPv 

Proof. Suppose w.l.o.g. that bo = 0, then Fq^^ is the convex hull of 
the origin in C^. and a finite set of points {vi, . . . ,vm} all of them 
satisfying {vi,v) = A^, < 0, hence lying in a hyperplane Hv not passing 
through the origin. 

The intersection of Fq^^ with a hyperplane parallel to Tiv will be 
nonempty if and only if A^, < g < 0. If this happens, the intersection 
will be the convex hull of 

{t-^i,--- ,7-^m}, 

Ay Ay 

which is equal to t"Qo„ = T~koPv □ 

Ai! ^ Ay 

In order to construct the Sylvester style matrix, we will take into ac- 
count whether the points lie in a translation of the primary cell or not. 
The first requirement we impose on 6 is that every point in {Q + 6)r\C 
must belong to the interior of a shifted maximal cell (primary or sec- 
ondary) . 

3.1.1. Points in the shifted primary cell. Proceed as in [CEl, CE2, 
Emil, Stu2]: choose generic lifting functions cDi, . . . , u;„, u, defined over 
Ai, . . . ,An, V{Q), in such a way that they produce a tight mixed co- 
herent decomposition of the Minkowski sum 

This implies that each n— dimensional cell JF in this decomposition 
equals 

Fi + F2 + . . . + F, + F, 
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where Fi is a cell in Aj,2i, -F is a cell in A^ and 

n = dim(Fi) + . . . + dim(F„) + dim(F). 

As a consequence, at least one of these dimensions is equal to 0. The 
row content of p G ({bo} + Qi + ... + Qn + Q + S) H C is a pair (i, a) 
defined as follows: 

• if p — 5 — bo lies in the cell J-' = Fi + . . . + Fn + F, and the set 
{j '■ 1 < J ^ n, dim(Fj) = 0} is not empty, let i be the largest 
index such that dim(Fj) = 0, and let Fi = {a}. 

• If dim(Fj) > 0, for all j, (because of the genericity of the lifting 
functions oJj, this implies that dim(F) = 0) then i := 0, and 
a := bo. 

We shall say that JF is mixed of type if the last item holds; otherwise, 
the cell shall be called non-mixed. 

Remark 3.4. The concepts of row content, mixed and non-mixed cells 
defined previously, appear with a slightly different meaning in [CEl, 
CE2, Stu2]. We shall discuss some relations between both definitions 
in example 3.3.3. 

We can now fill the rows of the matrix M indexed by those points 
p lying in ({bo} + Qi + ■ ■ ■ + Qn + Q + S) H £ as follows: for every 
p' G S, the entry indexed by {p,p') equals the coefficient of x^ in the 
expansion of the polynomial xP~"/j(x). Here, {i,a) is the row content 
of p. 

3.1.2. Points in the shifted secondary cells. Let v G M" such that 

(6) J-'v = -fo,D + Qlv + . . . + Qnv + Qv 

is a maximal cell in the polyhedral decomposition of Q. This implies 

that 

• dim(P^) = n- 1, Qi, = APj,, and Qi^ = kiP^, i = 0, 1, . . . ,n. 

• dim(Fo,i,) = n. 

Intersecting the shifted cell Fq^^ with hyperplanes parallel to P^ and 
using Lemma 3.3, it is easy to see that {J-'^ + 6) Ci C can be written as 
a disjoint union of sets of the type 

({\P^ + 6' + kiP^ + ... + knPv + \Pv + 6^nC, \e Q>o, 

which may be rearranged as follows: 

(7) {kiP^ + ... + knPv + XvPv + h) n £, A^ G Q>o, 
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where 6x := 6 + 6'. Here is where the recursion step comes: consider 
the v~ facet family 

(8) /i^(xi,... ,x„) = ^ Ci^ax" {i = l,... ,n). 

aeQi„nA 

Due to the fact that Qi^ = kiP^, has dimension exactly n — 1, it is 
straightforward to check that the family {Qi^ H ^j}j=i^...^„ is essential. 
This implies that the sparse resultant of the polynomials (8) is not 
constantly equal to one. We shall denote this sparse resultant as 

(9) Res^,(/i^,... ,/„J. 

In order to use the inductive hypothesis, we must decrease the dimen- 
sion of the supports with some care: let C^ C Che the lattice which is 
orthogonal to v, and denote by i3^i^+...+^„„ the afiine lattice generated 
by Ai^ H h Anv 

After a translation, we may suppose w.l.o.g. that G M" is a vertex 
of Py. This implies that £_4^^+...+^„^ is a sublattice of Cy, both having 
dimension n — 1, and we may consider the index [C^ : £_4^^+...+_4„^], 
which will be denoted by ind^. Let qi, ■ ■ ■ , qind^ be coset representatives 
for £^i^H_...+^„„ in Cy. 

For every p E C, there exists a unique j G {!,..., ind^} such that 
p G (qj + Cj[j^^+...+Anv) ® ^"^j so p may be written as pv + p" , the 
latter being an integer multiple of v, and p^ G q^ + £yii^+...+yi„^. Also, 
5a G £ ® Q = (£^,^+...+^„^ © Zw) (g) Q may be decomposed as 5^ + 5a„, 
where 81 (resp. 5xy) lies in Qv (resp. £^,^+...+^„^ (g)Q). 

If, in addition, p belongs to (7), we must have p" = 5^. This is due 
to the fact that kiPy + . . . + knPv + KPv + 8xy C -C^i^+.-.+y^^^ ® ^- So, 
p — p" = p — Sxy lies in 

{kiPy + ... + knPv + \Pv + hv) n (qj + >C^i„+...+yi„J . 

Finally, set (5j„ := 5\y — (\j G £^^^+...+^„^®Q. Now it is straightforward 
to check that a point p belongs to (7), if and only if, there exists 
j = 1, . . . , ind„, such that p := p — 5^ ~ 1i ^i^^ i'^ 

(10) (/ciP^ + . . . + A;„P^ + XyPy + (5^^) n £^i„+...+^„„. 

We have decreased dimension, and may use the inductive hypothesis 
in order to compute the resultant (9) using £_4^^+...+^^^ instead of C 
and 5jy instead of 5. It turns out that, for every A^ > 0, there will 
be ind^; square matrices of the form Ml^^x^ indexed by the points p 
lying in (10). Using the monomial bijection p = p + 5^ + q,-, we can 
relabel rows and columns of these matrices with the points of (7). The 
determinant of each Mj,^At, will be a non trivial multiple of (9), and each 
of these determinants will have the same degree as Res„(/i^, . . . , fn^) 
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in the coefficients of /i^ (/i^ will play here the role of /o in the previous 
step). 

Again, as in the primary cell, each of these matrices has in the row 
corresponding to a point p, the coordinates of the monomial expansion 
of x^~" fivi^) fo^ some a, i. In order to fill the row of M corresponding 
to p, we proceed as before: the entry indexed by {p,p') equals the 
coefficient of x^ in the expansion of the polynomial x^'"" fi{x). 

In order to finish the algorithm properly, the reader should check 
that, in the case n = 1, this procedure constructs a classical "Sylvester 
style matrix" ([Syl, Mac], see also Example 3.3.1) for two polynomials 
in one variable. If the supports generate the affine lattice Z, the matrix 
will be indexed by a set of monomials of the type 

{xl,x'l+\...,xl+'} aeZ,seN. 

Remark 3.5. Observe that, at each step of the recursion, we need to 
impose some conditions on the different 6^s in order to guarantee that 
all integer points are in the interior of a cell in each intermediate step. 
This happens for 6 generic. 

Remark 3.6. Let d^ be the v— lattice diameter of the cell JF^ which is 
defined as follows: 

(11) dv := max(m, f) — min(m, f). 

It is straightforward to check that, the number of matrices of the 
type Mj,^At, is exactly dv := d^mA^. Besides, due to the fact that 
Fi, F^i ■ ■ ■ iFni F are facets associated to v, the difference (11) may 
be actually computed as 

(12) max {m,v) — min {m,v). 

Remark 3.7. At each step of the recursion, the shifted secondary cells 
of the previous step are partitioned in such a way that their integer 
points are distributed into new primary and secondary cells. The new 
primary is again subdivided into mixed and unmixed cells. We shall 
keep track of this information, so a point will be said to be in a mixed 
cell of type i, if it belongs to a mixed cell which appeared at step i + 1. 
It is easy to see that, if a point is in a mixed cell of type i, then the row 
it indicates contains some coefficients of the expansion of a multiple of 

/.. 

It is also clear that, at the end of the recursion, each point in S has 
associated a row content. 
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3.2. Generalized Macaulay Formula. Now we are ready to state 
and prove the central result of this section. Before doing that, note 
that det(M) is well defined up to sign, because we have not given any 
order among the elements of £. Also, Res^(/o, . . . , /«) is well defined 
up to sign. So, the following statement will be true up to sign. 

Theorem 3.8. M is a generically nonsingular Sylvester style matrix. 
Moreover, we have the following formula "a la Macaulay" : 

det (M) = Res^(/o, ...,/„) det (E) , 

where E is the square submatrix of M formed by omiting all rows and 
columns indexed by points lying in mixed cells. 

Remark 3.9. It is easy to see that E does not contain coefficients of /o, 
so as an inmediate corollary of Theorem 3.8 we get that det (M) has 
the same degree as Res_4(/o, . . . , /„) in the coefficients of /q. Replacing 
the role played by /o with any /j, i = 1, . . . , n, we also have a formula 
for computing Res^(/o, . . . , /«) as the gcd of n + 1 determinants (see 
also [CEl, CE2, EM]). 

Proof. First of all, we will prove that, for a given point p G £, if p' 
belongs to the support of x^'"" fi, then p' must also be a point of S. 
Here, the pair (i, a) is the row content of the point p. This will imply 
that M is a Sylvester style matrix. 

Two different scenarios must be considered: 

• If the point belongs to the shifted primary cell, proceeding as in 
[CEl, CE2], it is easy to see that this happens. 

• If the point belongs to a shifted secondary cell, let us say ^F^, 
because of (6), x^'"" fi^ has its support contained in 

This, combined with the fact that, in secondary cells, i is always 
bigger than 0, implies that the support of x^~"'fi is contained in 

Fo,v + Qlv + • • • + Qi-lv + Qi + Qi+lv + • • • Qnv ~\~ Qv + ^i 

the latter set being a subset oi Q + 5. From here, the claim follows 
straightforwardly. 

Once we know that M is a Sylvester style matrix, it is easy to see 
that Res^(/o, . . . , /«) divides det (M) using the standard argumenta- 
tion given in [CEl, CE2, Stu2]. 

We shall regard det (M) as a polynomial in 'L[ci^a\ and will prove that 
it is not identically zero by showing that its highest term with respect 
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to some monomial order is nonzero. Explicitly, we shall prove that, for 
the vector 

(13) uJ:={uJo,uJi,... ,uJn)eW^, 
where the Ui were defined in (4), 

(14) initzj {det (M)) y^O. 

It is easy to see that, writing det (M) as a polynomial in co,bo with 
coefficients in Z[cia \ {cobo}]) the leading term of this polynomial is 
(14). 

Due to the special role that /o has played in the construction of the 
matrix, the number of integer points lying in shifted 0— mixed cells is 
equal to Mq := A4V{Qi, . . . , Qn), i-e. the degree of Res^(/o, . . . , /„) in 
the coefficients of /o ([HS, CEl, CE2, Stu2]). So, it is straightforward 
to check that 

d^gcoeff(/o) (Res^(/o' ■■■ Jn)) = degcoeff(/o) (det(M)) . 

This, combined with the special way in which we have lifted the poly- 
topes (i.e. just hfting the point bo), implies that 

det(M) _ initjj{det{M)) 

Res^(/o, ... , /„) initjj (Res^(/o, . . . , /„)) 

coefficient ofcj^^^in det(M) 



coefficient of cj^^^ in Res^(/o, ... , /„) 

In order to prove the theorem, we shall proceed as in [Mac], by showing 
that the numerator of this fraction is nonzero, and that the extraneous 
factor, i.e. the ratio, is det(E). The proof will be again by recurrence 
on n. 

The basic case (n = 1) is completely contained in the classical for- 
mulas given by Macaulay in [Mac]. In the general case, in order to 
compute init jj {det (M)) , we proceed as in [CEl, CE2, Stu2]: replace 
the polynomials (2) by the following deformed family: 

(15) fi,^ := J2 c^.'^ t'^'^'^x", i = 0,l,...,n, 

aeAi 

and consider the deformed matrix M (cj^a^'^' *■"•*) • 
Remark 3.10. Actually, 

fi,uj = fi i>l- 

It is easy to see that initu {det (M)) is the leading coefficient of the 
determinant of the matrix M(cj^a^'^ ) regarded as a polynomial in t. 
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For every p & S, let h{p) be the biggest rational number such that 
(p - 5, h{p)) eQn = Qo,n + Qi,n + ■■■ + Qn,n + Qn- 
The following observations will be useful later: 

Lemma 3.11. The function h verifies 

1. < h{p) < 1, for all pes. 

2. h{p) = 1 if and only if p ~ 5 lies in the primary cell. 

3. If p — 6 and q — 6 are in the same secondary cell, say JF^, then 

h{p) = h{q) ^^ {p^v) = {q,v). 

4. If p E {J-'v + S) n C, the row content of p is the pair {i,a), and 
v' j^ fiv, /i > 0, then 

{{p-5- a, h{p)) + Qi^n) n Qn(.',r) = 0, VZ' G R<o; 

here, Qn(v',i') C W^^^ is the face of Qq determined by {v',l'). 

Proof of the Lemma: The first two statements are obvious. The third 
assertion holds straightforwardly just noting that JF^, is the projection 
of Qn(v,i)j where / is a negative real number, and 

p - 5, q - 6 e Tv ^^ (P - S, h{p)) ,{q~6, h{q)) G Qn{v,i), 

hence 

{ip-5,hip)),ivj)) = {iq-6,hiq)),iv,l)). 

Lastly, let / G M<o such that Qn{v,i) projects bijectively onto J-'^. Due 
to the generic conditions imposed on 6, it turns out that the point 
Pn '■= {p — S, h{p)) belongs to the relative interior of Qq.(v,1)- If v' = fiv, 
and I' = fil, with /i < 0, this would imply /' G ]R>o, which is not of 
interest for us, so we can suppose w.l.o.g. that (f ', /') is not parallel to 
(f , /). Then, one can slightly displace the point pn inside Qn{v,i), in the 
direction of the orthogonal projection of —{v',l') over the hyperplane 
{{x, {v,l)) = 0} C R"+^. After this displacement, all the points in the 
shifted Qi will still lie in Q^. This is due to the convexity argument 
given in [CEl, CE2], which states that for every point q^ lying in 

qn - (a, 0) + Qi^n C Qq. 
So, if the statement of the Lemma were not true, the common points 
in the intersection with Qn{v',i') would not belong to Qn after the dis- 
placement, which is impossible. 

D 

We shall use the convexity argument given in [CEl, CE2, Stu2] as 
follows: for each p E S, we multiply every element in the row indexed 
by p by t^Pi^'^iw ^ where the row content of p is (i,a). Let us call the 
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matrix obtained in this way M'(t). It is easy to see that the leading 
coefficient of det(M'(t)) (as a polynomial in t) is initjj{det{M.)) . 

Proposition 3.12. Let < 71 < 72 < . . . < 7Ar = 1 be the differ- 
ent values for h{p) as p ranges in £. Then, the leading coefficient of 
det(M'(t)) (as a polynomial in t) factorizes as 



N 

(16) JJdet 

j=i 



'^3) 1 



where M^ is the square submatrix of M made by choosing all rows and 
columns indexed by points p such that h{p) = 7j. 

Moreover, the product (16) can be reorganised as follows: 

(17) det(M;v) n (det(Mi) . . . det(M^)) , 

where the product is taken over all secondary cells JF^ and M^ is one 
of the matrices Mj,^At, defined in Section 3.1.2, i.e. a Sylvester style 
matrix for the sparse resultant associated to the family (8). 

Proof of the Proposition: We will compute the leading term of det(M'(t)) 
by searching, in each column, the highest power of t appearing in that 
column, replacing by all the entries which have not this highest power, 
and computing the determinant of the modified matrix. We shall call 
this matrix mod{Wl'{t)). 

It is straightforward to check that, in M'(t), the highest power ap- 
pearing in the column indexed by p' is exactly h{p'). So, in mod(M' (t)) , 
in the columns indexed by those p' such that h{p') = 71 there cannot 
be nonzero entries in the rows indexed by those p such that h{p) 7^ 71; 
otherwise h{p) > 71 and this will imply that the power of t appearing 
in the {p,p') place of the matrix will be strictly bigger than 71, which 
is impossible. 

This implies that mod{Wl'{t)), after division by some power of t, 
and ordering its rows and columns by putting the points p' such that 
h{p') = 7i at the beginning, has the following structure: 



Ml A 
B 
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Repeating the argument recursively, it turns out that the structure of 
mo(i(M'(t)) is triangular as follows: 



/Ml * 
M2 



V 



* 



UN 



J 



and the first part of the proposition holds straightforwardly. 

In order to prove (17), we fix j < A^. Then, the points p such that 
h{p) = 7j < 1, lie in shifted secondary cells, let's say jFj,^ +6, . . . J-'vj^+6, 
and one can arrange the rows and columns of the matrix Mj such that 
all the points in JF^,^ + S appear at the beginning, the points in jF^j + ^ 
immediately after, and so on. 

First of all, we will show that 



/ 
















\ 



where 



Hj,vk) 



V . 

denotes the submatrix of 



whose rows and columns 



are indexed by points in the cell JF^^ + 5. 

In order to do this, consider the deformed family (15). It is straight- 
forward to check that 

supp (t^(P)-'^>(")xP-"/.,.) = {p-a, h{p) - u.ia)) + Q,,n C Qn + {S, 0) 

where, as usual, {i,a) is the row content of p. Moreover, the point 

{p - a, h{p) - uJi{a)) = {p-a, h{p)) 

belongs to the facet of the shifted polytope Qn + {S, 0) determined by 
and inward normal vector of the type {vk, I), I & IR<o- 

Because of the last item of Lemma 3.11, there cannot be nonzero 
coefficients corresponding to the expansion of the polynomial 

(19) t^(p)~^'(«)xP-"/*,.. 

whose multidegree in {x, t) lies on the boundary of Qn + (5, 0) other 
than those in the facet determined by {vk,l)- This implies that, if the 
point {q, s) is an exponent arising in the expansion of (19) and q does 
not belong to the shifted secondary cell J-'y^ + 6, then (g, s) must be 
an interior point of Qn + (5,0), and this implies that s < h{q). Due 
to the remark made at the beginning of this proof, it turns out that 
the element indexed by {p,q) in mod(M'{t)) is zero. This gives the 
structure stated of M^. 
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Now, recalling Remark 3.6, and using the third item of Lemma 3.11, 
it turns out that M(j^^^) is actually a matrix of type M^^^x^ . Moreover, 
all matrices of this type appear in this way, so the proposition holds 
straightforwardly. D 

Let us return to the proof of the Theorem. By the inductive hypoth- 
esis, det(M^) 7^ 0, \/v, i, and using the special lifting Coi, . . . u;„, u in the 
primary cell (see Subsection 3.1.1), it is easy to see that det(M7v) 7^ 0. 
Moreover, 

(20) det(M^) = c^^^ det(E^) ^ 0, 

where Kn is the submatrix of Mat formed by all the rows and columns 
indexed by points in non- mixed cells (see [CEl, CE2]). This proves 
that det(M) ^ 0. 

Again by the inductive hypothesis, det(M^) equals Res„(/ij,, . . . , /„„) 
times a minor det(E^) made by choosing all rows and columns non- 
mixed in M^. This implies that 

(21) 

dv 

zmt^detiM) = c^l detiE^) J] (Res(/i„ . . . , U,)f- J] J] det(E^). 

By selecting all rows and columns indexed by non-mixed points in M, 
a modified version of Proposition 3.12 holds for the matrix E, which 
also has a block structure, and we get 

(22) det(E) = det(E^) JJ JJ det(E;). 

The proof of the theorem will be complete if we show that 

^mt^(Res^(/o, ...,/„))= <^„ J] (I^es„(/i„, . . . , /„J)^- . 

Using Theorem 4.1 in [Stu2], we have that 

imtjj(Res^(/o, ... ,/„)) = ± JJ (Res^(/o|Fo, • • • , /n|F„))'^^ , 

where JF runs over all facets of the coherent mixed decomposition given 
by cJ in the Minkowski sum 

Q:=Qo + Qi + ... + Qn- 
Explicitly we have that 

^=Fo + ...+Fn 
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with Fi = conv{A'i), A'^ C Ai, 

Ji\Fi -^ / ^ Ci^aX , 

Resjp is the sparse resuhant associated to the data {A'q, . . . ,An), and 
d^ equals the unique integer such that (ReSjp-(/o|Fo, • • • ;/n|F„)) ^ has 
total degree Zi=o MV{Fo, . . . , Fj_i, F,+i, . . . , F„). 

It is easy to see that the coherent mixed decomposition induced by 
Co over Q is similar to the one induced over Q. More precisely, we get 
a big primary cell of the form ho + Qi + ... + Qn, and for each v such 
that J-'y is a secondary cell of Q, there is a secondary cell JF„ in Q. 
Moreover, the same analysis made in 3.1 says that these are all the 
cells of maximal dimension in Q. 

For the primary cell, we get that 

flip, = fi i = l,... ,n. 

This implies that the unique essential set in the data (^q, . . . ,A'n) is 
just the singleton {bo}, so Res^(/o|Fo, • • • , UIfJ = Co,bo and d^ = 
MV{Q,,... ,g„) = Mo. 

Now, take a vector v such that JF^, is a cell of maximal dimension. 
This implies that /o|fo has an n-dimensional support and 

fi\F, = fiv i = l,...,n. 

This implies that, in this case, the unique essential set is {AinQi^}i=i^,„ ,„ 
and Res^J/olFo, • • • , IuIfJ = Res^(/i^, . . . , /„J. 

It remains to prove that d-p = dy. To begin with, observe that 

MV{F,, . . . , F„) = MV{Qu, . . . , Q„J = 

due to the fact that the supports lie in a hyperplane, so the Minkowski 
sum of any subfamily of {Fi, . . . , F„} does not have positive n-volume. 
In order to compute the other numbers involved in the computation 
of d^ , we shall use the recursive relation satisfied by the mixed volume 
([Ber" CLO]): 

MVn {Fo, Fi, . . . , Fl_^, F,+i, . . . , F„) = 
Z., ap.iv') MV'^_, m),,, (F,_i),,, (F^+i),,, . . . , (F„),0 , 

the summation being taken over all v' such that P^' is a facet. Here, 
-^^n-i ii'^i)v') denotes the normalized mixed volume with respect to 
the hyperplane v' G C orthogonal to v', and 

apoiv') := — min {m,v'). 

m£Fo 
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Using the fact that Fj = Qi^, it turns out that 

MV'^_, ((Fi),,, (F,_i),,, (F,+i),,, . . . , (F„)„,) = 0, 
unless v' = V OT v' = —v. Hence, we have that 

MVn{Fo,F,,... ,Fi_,,Fi,... ,F„) = 
(apM + aF,i-v)) -MV^Li ((^i)., iFi^i)v, (i^m). • • • , (i^n).) 
and therefore, due to the fact that 



deg(Res,) = T^Y1 -^^"-1 ((^i)- (^'-i)- (^'+i)- • • • ' (^")-)' 

we get that 

d-^^ = {apoiv) + o-Foi-v)) indt, = dv 
as claimed. D 

Corollary 3.13. For every i = 0,1, . . . ,n, the number of points lying 
in a shifted mixed cell of type i is exactly A4V{Qo, . . . , Qi-i, Qi+i, ■ ■ ■ , Qn) 

Proof. It is straightforward to check the following equalities: 

degpQgg . (det(E)) = #{non-mixed points of type i}, 
degpQgg „ (det(M)) = #{non-mixed points of type i} 

+#{mixed points of type i}. 

U 

Remark 3.14. We required bo to be a vertex of ^o just in order to 
decrease the number of secondary cells in the algorithm, but this con- 
dition is not used in the proof. As an easy consequence, we get that 

// the family {/i}j=o,...,n is essential then, every generic coefficient 
Ci,a of the input support appears in Res^(/o, . . . , /„) with highest power 
MV{Qo, ... , Qi-i, Qi+i, ... ,Qn)- 

Remark 3.15. We have used the lifting algorithm of Canny and Emiris 
in primary cells in order to break ties, but the Theorem holds just 
provided that one can construct recursively Sylvester style matrices 
having non-zero determinant and the same degree as the resultant in 
the coefficients of any /j. This has been already noted by Macaulay in 
the classical case (see [Mac, Section 6a]). 
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3.3. Examples. 



3.3.1. The Onedimensional Case. Set 

A := {0,2,4}, A:={4,8}. 

Here, the affine lattice C equals 2 Z, and the polytope P is the unit seg- 
ment [0, 1]. Set A := |, (5 := |, and bo := 0. Then, it is straightforward 
to check that: 

. £ = [4 + i, 14 + |] n 2 Z = {6, 8, 10, 12, 14}. 

• The points lying in the shifted primary cell are 6, 8 and 10; the 
other points belong to the (unique) shifted secondary cell. 

• There is a unique rule for filling the rows of the matrix M cor- 
responding to the points lying in the shifted secondary cell. We 
have that 



x'^ /i (mixed) 
x^ /i (mixed) . 



Although there are infinitely many different lifting functions Coi, Co 
over Ai and V{Q) = {0, |}, they cannot produce more than two 
different tight coherent mixed decomposition of the segment 



[4,8 + ^] 



Qi + Q. 



Explicitly, we get the following cases: 
1. 



x^ 


1 — > 


x^ 


/o 


(mixed) 


x« 


1 — ¥ 


x« 


/o 


(mixed) 


^lU 


1 — ¥ 


x' 


/l 


(non- mixed). 



which corresponds to liftings which give the same partition 
than 



c^i(4,8 

'^(0,1; 



(0,1) 
(0,0); 



2. 




corresponding to 



c:'i(4,8 

^(0,1; 



(0,0) 
(0,1). 
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Observe that both cases give essentially the same matrix, the only 
difference is that the rows are indexed differently. Setting 



/o 
/i 



(23) 

we get the following matrix 



ex 



a + hx\ 



15 



(ah 

a 

d 

rf 
y 



\ 

c 

e 
del 



The matrix E here consists by the element e which appears in the third 
row and column of M. Expanding the determinant by the last column, 
we have that 

/ a 6 c \ 

a 6 c 

rf e 

y rf e y 

and we can easily see that the determinant of matrix on the right hand 
side corresponds to the Sylvester resultant for the bivariate family given 
by (23). 



detf 



Remark 3.16. It is not hard to see that, for every bivariate family, 
every A, bo and generic 5, we get a scenario similar to this example. 
More precisely, the algorithm given in the previous section produces 
the same Sylvester style matrix given by Macaulay in [Mac]. 

3.3.2. Let us compute the resultant of the Example 2.4. Take A = 1, so 
Q will be the unit square. Set also bo = (0, 0). In this case, C = iP, the 
primary cell equals to 3 times [0, 1] x [0, 1] and there are two secondary 
cells, corresponding to the inward vectors (—1,0) and (0, —1). 

Setting S = (|, |), we get that £ consists of sixteen points, nine of 
them lying in the shifted primary cell. Explicitly, we have that 

£ = {(a,b) el? : 1 <a,6<4}. 

In this case, the facet polynomials {/ij,,/2t,} can be regarded as two 
polynomials in one variable. So, the rows indexed by points lying in 
shifted secondary cells may be filled in the classical Sylvester style as 
in the previous example. In order to fill the rows corresponding to 
points lying in the shifted primary cell, we will use the following lifting 
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functions on the ordered sets At = V{Q) = {(0, 0), (1, 0) (0, 1) (1, 1)} 



a) 



(0,1,1,2); 

(0,0,7,7); 

(0,14,0,14). 



Observe that the hfting functions are actually the restriction of a linear 
function on R^. Using any algorithm for computing convex hulls like 
the one given by Emiris in [Emi2], we obtain the following description 

for M: 



row 


coefficients of 




cell 


type 


0:10:2 


X1X2 /i 




primary 


non-mixed 


xlx2 


xlx2 /o 




primary 


mixed 


X\X2 


xlx2 /l 




primary 


non-mixed 


xixl 


3^1X2 /2 




primary 


non-mixed 


2 2 

X-i Xn 


XIX2 /2 




primary 


non-mixed 




2^l2^2 /o 




primary 


mixed 


X\X2 


Xixl /2 




primary 


non-mixed 




a^ia^i /2 




primary 


non-mixed 




■^l'^2 /l 




primary 


non-mixed 




■^l'^2 /l 


(-1 


0) — secondary 


mixed 


12 


a;?a;i /a 


(-1 


0) — secondary 


non-mixed 


X\X2 


a;?X2 /2 


(-1 


0) — secondary 


mixed 


4 
X\X2 


xixl /2 


(0,- 


-1) — secondary 


mixed 


12 


x^xi /a 


(0,- 


-1) — secondary 


non-mixed 


12 


a;?a;i /a 


(0,- 


-1) — secondary 


non-mixed 


12 


X\X2 /l 


(0,- 


-1) — secondary 


mixed 
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Writing fi = ai + biXi + CiX2 + diXiX2 and ordering the monomials as 
in the table, we get that 



/ Ol 


bi 





Cl 


dl 
































\ 





ao 


^0 





Co 


do 



































ai 


bi 





Cl 


di 
































02 


b2 





C2 


d2 






































02 


62 





C2 


d2 















































oo 








Co 


da 


bo 


























02 


62 





C2 


d2 






































02 


62 





C2 


d2 



































01 


bi 





Cl 


di 






































01 








Cl 


di 


bi 
































02 








C2 


d2 


b2 























«2 








C2 














d2 


b2 
































02 


b2 














C2 


d2 





























02 


62 














C2 


d2 





























02 


b2 














C2 


d2 


V 























01 


bi 














Cl 


di 1 



and 



E 



/aiOcirfiO 0\ 






hi 





Cl 




















a2 





C2 


d2 























62 





C2 


























02 


62 


C2 


d2 























02 





C2 


d2 




















fll 





Cl 


d. 





























C2 


62 























02 


&2 





C2 


^2 



V 



02 C2 / 



With the aid of Maple, we can check that 



det( 



±Res^(/o,/i,/„)det(E) 



and that det(E) factorizes as 



-C2(-cia2 + aiC2)b2{cid2 - diC2){-h2Ci + 61C2) 
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In this easy example, one can also check that the leading term of det(M) 
as a polynomial in ao is the determinant of the following matrix: 

(24) 



2i 


bi 





Cl 


d, 
































\ 





ao 















































ai 


h 





Cl 


d. 
































12 


&2 





C2 


d2 






































02 


&2 





C2 


d2 















































oo 









































02 


&2 





C2 


d2 






































02 


&2 





C2 


d2 



































Ol 


h 





Cl 


d. 


















































d. 


bi 












































d2 


b2 















































d2 


&2 


















































C2 


d2 















































C2 


d2 















































C2 


d2 












































Cl 


di ) 



and we can recognize in this matrix the block structure stated in Propo- 
sition 3.12. Explicitly, we have a 9 x 9 big block coming from the pri- 
mary cell, and two blocks arising from the secondary cells, of size 3x3 
and 4x4 respectively. Computing the determinant of (24), we get: 

-aoC2(-cia2 + aiC2)b2{-C2di + d2Cif{bid2 - &2'^i)(c2&i - C162), 
and we can check that 



and 



can check that 

det(E) = cl{-cia2 + aiC2)&2(-C2C?i + c?2Ci)(c2fei - C162), 

Res(o-i)(ciX2 + (iia;iX2,C2a;2 + d2XiX2) = {-C2di + d2Ci), 
Kes(^ifl){biXi + diXiX2, b2X2 + d2XiX2) = {bid2 - &2'^i) 



as expected. 



3.3.3. We want to compute the sparse resultant of the family 

/o = «1 + Ci2Xl + 03X2, 

/i = bi + b2Xi + 63X2 + bixl + 65X1X2 + bexl, 

/2 = Cl + C2X1 + C3X2 + C4X1 + C5X1X2 + CqxI + 
+0^x1 + Csxlx2 + CgXiXJ + CiqxI. 

Here, Ai = {{a,b) eNl : a + b <i + l}, i = 0,1,2, the lattice S coin- 
cides with Z^ and the polytopes Qi are integer multiples of the standard 



(25) 
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simplex 5*2. The sparse resultant coincides with the classical resultant 
of three homogeneous polynomials of degrees 1,2 and 3 respectively, 
whose afiinizations are the fi (see [CLO, Mac]). 

In order to compute this resultant, we set Q := and bo := (0,0). 
The primary cell equals 5 5*2, and there is a unique secondary cell, 
corresponding to the vector (—1, —1). 

Setting 6 := (e, e), with 1 :^ e > 0, we have that S has 15 monomials, 
ten of them lying in the shifted primary cell. We use the following lifting 
function, defined on the vertices of each Qi, i = 1,2 ordered as follows: 
{(0,0), {k,0), (0, A;)} where k = 2,3, and extended to the rest of the 
points of the input support by linearity: 



UJi 

C02 



(1,0,1) 

:i,i,o). 



In this case, the subdivision does not depend on the value of u, the 
lifting function over the unique vertex of Q. As in the previous ex- 
amples, the secondary cell will be filled in such a way that the facet 
resultants will be computed using the classical Sylvester formula for 
bivariate polynomials. 
Explicitly, we get: 



row 


coefficients of 




cell 


type 


JbAJbiy 


XxX2 f\ 




primary 


non-mixed 


xix\ 


^\^2 /l 




primary 


non-mixed 


xjxl 


X\X2 /l 




primary 


non-mixed 


X\X2 


3^1X2 /2 




primary 


non-mixed 


X\X2 


X\X2 /2 


^ -*-? 


— 1) — secondary 


mixed 


xxx\ 


X\X2 f\ 


' — 1, 


— 1) — secondary 


mixed 


x\x\ 


^V^2 /l 


^ -*-; 


— 1) — secondary 


mixed 


x\x\ 


X\X2 /l 


^ -*-? 


— 1) — secondary 


mixed 


x\x\ 


Xxx\ /2 


^ ^^ 


— 1) — secondary 


mixed 


X\X2 


X\X2 /o 




primary 


mixed 


Uj\Ujn 


xix\ /o 




primary 


mixed 


x\x2 


A^2 /o 




primary 


mixed 


x\x2 


X\^2 /o 




primary 


mixed 


JLa Jb'^ 


x\x\ /o 




primary 


mixed 


JL -\ JL^y 


^\^2 Jo 




primary 


mixed 
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and, indexing the matrix M using this order, we get 



/ h 


























61 


^3 


b2 


65 


h 


\ 


bs 


h 


h 























bi 





b2 





64 








h 


64 























61 


bs 


b2 


b5 


C6 


ClO 


eg 


C7 

















Cl 


C3 


C2 


C5 


C4 


C8 








C6 


C4 


C7 





ClO 


eg 


C8 








Cl 


C3 


C2 


C5 


h 


^3 


b2 








^6 


b5 


64 





























b3 











ba 


b5 


bA 











bi 





62 











62 


64 








ba 


b5 














bi 


63 


C3 


C6 


C5 








ClO 


eg 


C8 


cr 





Cl 





C2 





C4 





























01 


03 


02 











03 





























01 





02 









































ai 


03 


02 














02 





























01 


03 























03 


02 

















01 


\ 





03 





























ai 





02 / 



In this case, the extraneous factor E is the submatrix made by choosing 
the first four rows and columns of M, and its determinant equals 

be (c7 bl - 64 be cg + 64 65 Cio). 

Observe that all points lying in shifted non-mixed cells are actually in 
the shifted the primary cell, so there is a priori no significant block 
structure in this matrix. 

Curiosity: The matrix constructed here is exactly the same given by 
Canny and Emiris in the last section of [CE2] (see also section 3.1.4 
of [Emil]) in order to show that, in their construction, the extraneous 
factor is not always the determinant of the minor formed by choosing 
all rows and columns indexed by non-mixed points. Of course, they 
work with another definition of non-mixed cells! 

To be more precise, their algorithm is not recursive. They produce 
this matrix by applying a lifting algorithm to the polytopes Qo, Qi and 
Q2 using the linear functions /q := iO'^Xi + 10'^X2, h := lO^Xi, I2 : = 
10^Xi-|-10a;2. By taking the lower hull of the lifted polytopes {Qi, kiQi)), 
a tight coherent mixed decomposition of the Newton polytope Q0+Q1 + 
Q2 holds. Using 6 as before, the same matrix is constructed but the 
points lying in the shifted non-mixed cells are in correspondence with 
the monomials 

iJyJ^JyO. JulJbn^ JjAjb<2i Jb'\Jb'2 \ • 

By taking the determinant of the submatrix of M made choosing the 
rows and columns indexed by these monomials, we get b\c-j. 
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3.4. An Overview of Macaulay's Classical Formulas. In this sec- 
tion, we will see how the formulas given by Macaulay in [Mac] can 
be recovered with our methods. In order to have a notation similar 
to Macaulay's original paper, we shall deal with n generic polynomi- 
als in n — 1 variables Xi^ . . . ,x„_i of total degree less or equal than 
mi,... ,m„ respectively. In our terminology, the input supports are 
integer multiples of the standard simplex S'„_i. More precisely, the 
polynomials will be denoted as Ci, C2, • • • , Cn and 

supp{Ci) =miSn-i. 

Actually, Macaulay worked with n homogeneous polynomials in n vari- 
ables, but it turns out that the homogeneous resultant of these poly- 
nomials coincides with the sparse resultant associated to the supports 
Ci as in Section 2. 

Remark 3.17. Setting t„ := Yll=ih^i ~ l)? ^^"^ using the algorithm 
given by Canny and Emiris, it is also possible to recover Macaulay's 
classical formula in degree t = t„ + 1. Moreover, the extraneous factor 
in Macaulay's original formulation is exactly the minor formed by using 
all rows and columns indexed by points lying in shifted non-mixed cells 
which they get with their methods (see [CE2] and [Emil]). 

In section 3 of [Mac], Macaulay constructed a Sylvester style matrix 
whose determinant is denoted by D{n, t), for every t G N>o. The matrix 
has its rows and columns indexed by all monomials of total degree less 
or equal than t. For every t > t„, it turns out that D{n, t) is a nonzero 
multiple of Res^(Ci, . . . , C„) denoted in that paper as R{n^ t). ^ 

His construction is recursive in the following sense: let v be the 
inward normal vector (— 1, — 1, . . . , — 1); it is easy to see that 



(26) Ci,^, C2ut • • • 5 C*9 



n— It 



is a family of n— 1 homogeneous polynomials in the variables Xi, . . . , x„_i, 
so after setting Xn~i equal to one, it may be regarded as a sparse family 
with support in 

f^l Sn-2, 1^2 Sn-2i • • • 5 ^^n-1 Sn~2i 

SO D{n — 1, j) will be the determinant of a Sylvester matrix made with 
the same rules as D{n^ j) but with n—2 variables, using the polynomials 
(26). 



^For lower values of t, Macaulay also proposed a matrix with similar properties, 
but not of Sylvester type. 
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In section 5 of [Mac], he established the following theorem: 

^ ' R{n,t) 11 R{n-l,t-j) 11 ^ ^' 

which is the cornerstone of the main result given in section 6, namely, 
that R{n, t) can be recovered as the quotient of D{n, t) by the minor of 
obtained by omitting all rows and columns corresponding to monomials 
reduced in all the variables. 

Here, a monomial x" = x"^ . . . , x'^Ti is said to be reduced if there 

exists a unique i G {1, . . . , n} such that x™' divides x" xiT . 

Macaulay gave also a recursive structure of the extraneous factor (in 
his notation, A{n,t)) as follows: 

(28) A(n, t)= H A{n-l,t- j) J] D{n - 1, k) 

j=0 k=l 

(see section 6 of [Mac]). 

In order to see this construction in light of the results presented in 
the previous section, extra care must be taken, because Macaulay's 
construction produces a determinant which has the same degree as the 
resultant in the coefficients of C„, so we shall modify our algorithm in 
such a way that the role of /o is played by C„; similarly, the role of /i 
will be played by C„_i, and so on. 

The polytope P will be the standard simplex Sn-i- Given t > tn, in 
order to define Q, set 

X := t — tn — 1, 
so Q will be equal to 

(mi+. . .+m„+A)S'„_i = {(mi, . . . , m„__i) G M"""^ : < m^ < t+n-1}, 

and C = Z"^^. By taking 6 = (e, . . . , e), with 1 ^ e > 0, we get that 

^ = {(«!,... ,a„_i) gZ"-^ : l<a„ ^a, <t + n-l}. 

i 

Setting 

.29) ^ ^ {/3eN^l: EA<t} 

(«!,... ,a„_i) h^ (ai - 1,... ,a„_i - 1), 

we get a bijection between our support and the one used by Macaulay 
for computing the matrix whose determinant is D{n,t). By choosing 
bo = (0, . . . ,0), we can check that there is only one secondary cell, 
associated to the vector u, whose z/— diameter is exactly m„. 

Using lifting functions cJi, . . . , Un-i, Co as in section 8 of [CE2], it is 
possible to get a subdivision of the primary cell such that the points 
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lying in the 0— mixed cells are bijectively associated (with the bijection 
given in (29)) with those reduced monomials which are divisible by 
x™". Lifting functions with the same properties should be used in the 
recursive steps. 

Now, comparing equations (28) and (22), it is not hard to check that 

m„-l 

n A(n-i,t-j)=nndet(n)- 

i=0 T^ i 



Hence, 



t—nin 



det(E^) = JJ D{n-l,k). 



k=l 



4. The general case 

We will extend here the results of the previous section by assuming 
only that the family {Ai}o<i<n is essential, without any other con- 
dition on the supports. Observe that this hypothesis ensures that 
Res^(/o, • • • , /n) is non trivial. 

As before, set 

Qi = conv{Ai), i = 0, . . . , n, 

and consider the Minkowski sum Qo+. . ■+Qn- As the family of supports 
is essential, this polytope must be n-dimensional. So, £ir = M". 

In order to give more generality to our algorithm, let us consider 
a "generic" polytope Q (in a sense which will be determined later, 
see for instance Remark 4.4), with vertices in £ Q. The algorithm 
will produce a Sylvester style matrix whose rows and columns will be 
indexed by the elements of 

(30) £:={Qo + Qi + ... + Qn + Q)r]C, 

and whose determinant will be a non zero multiple of the sparse resul- 
tant. The extraneous factor will be a minor of this matrix. 

As before, the algorithm will be recursive on the dimension of the 
polytope 

Q:=Qo + Qi + --- + Qn + Q, 

and in its intermediate steps will use the subdivision technique of 
Canny and Emiris (see [CEl, CE2, Stu2]). 
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4.1. Construction of the Sylvester matrix. Let Q be polytope as 
before, V{Q) C £ (g> Q be the set of vertices of Q. 

Let us pick a vertex bo G ^O) and consider the same lifting functions 
oJi, uj, defined in (4). Consider also 

Qi,n ■= conv{{a,uJi{a)) : a e Ai}, 
Qn := conv{{b,ij{b)) : beV{Q)}, 

and the coherent mixed decomposition Aj ^^ (resp. Aj^) of the polytopes 
Qi (resp. Q) given by projecting the upper envelope of Qi^Q (resp. Qn). 
In the decomposition of Q, each cell is of the form 

where Fj (resp. F) is a cell in Aj^^, (resp. A^). 

Our lifting functions are again mostly trivial. So, we can compute 
the cells which will appear in the subdivision: in Ai^q, . . . ,An,n,An 
we may have two types of cells as in the previous section. 

The decomposition Ao,q has the following cells: the face determined 
by (0, —1) is just the singleton {bo}, and it is straightforward to check 
that, for every k > 0, every /c-dimensional cell in the decomposition is 
a A;-dimensional face of Qo which contains {bo}, or the convex hull of 
bo and a (A; — l)-dimensional face of Qq which does not contain this 
point. We shall denote it by Fq^^, where v is some interior primitive 
normal vector of that face. 

So, all cells in the decomposition of Q are as follows: 

• the primary cell: {bo} + Qi + . . . + Qn + Q- Observe that it has 
always dimension n due to the fact that the family {^i}i=o,i,--,n 
is essential. 

• -fo,D + Qiv + • • • + Qnv + Qv for some v G M". Those being n- 
dimensional will be called secondary, and they will be uniquely 
determined by their primitive inward vector v G M". 

As in the previous section, the Sylvester matrix will take into account 
if the points lie in the primary cell or not. The first assumption we will 
make on Q is that all points in Q fl £ must belong to the interior of a 
maximal cell (primary or secondary). 

4.LL Points in the primary cell. Proceed exactly as in 3. LI: choose 
generic hfting functions Ui, . . . ,Un,oj over en Ai, ... ,An, V{Q) re- 
spectively, in such a way that they produce a tight coherent mixed 
decomposition of the Minkowski sum 

Ql + ... + Qn + Q 
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by taking the upper envelope. Every n-dimensional cell T in the de- 
composition may be obtained as 

Fi + F2 + . . . + F„ + F, 

where Fj is a cell of l^i^;:^, F is a cell of A^, 

n = dim(Fi) + . . . + dim(F„) + dim(F) 

and at least one of these dimensions is zero. 

The concepts of row content, mixed cells of type 0, and the entries 
of M whose row coordinates are indexed by 

p e ({bo} + Qi + ... + Qn + Q)nC, 

are defined in the same way as in 3.1.1. 

4.1.2. Points in the secondary cell. Here, we cannot use the recursive 
step given in 3.1.2, so we must proceed with some care. Let v G M" 
the primitive inward normal vector of the n-dimensional cell 

(31) J-'v = Fo,v + Qlv + • • • + Qnv + Qv- 

Consider the family of polynomials associated to the f-facet 

(32) fi^{xi,... ,Xn) = ^ Q,aa;" (i = 1,... ,n). 

We cannot claim now that the family 

(33) {Q^,r]A^}^<i<n 

is essential, but certainly there always exists a subfamily of indices 
I < ii < i2 < ■ ■ ■ < ik ^ n, such that 

(34) Wv.nA,} 

is essential. It is not always true that (34) is the unique essential sub- 
family of (33), but we can assume w.l.o.g. that there exists an essential 
subfamily which is minimal with respect to the inclusion. Moreover, 
we may suppose that this subfamily is indexed as follows: 

(35) {Q^^nA^}l<^<k. 

This means that the sparse resultant of the polynomials /i^,, /2i,, ■ ■ ■ , fkv 
with respect to the supports (35) is non trivial. We shall denote it as 

(36) Res^(/i„,... ,/fc„). 

In order to mimic the inductive step of the previous section, suppose 
w.l.o.g. that G Aiy, i = 1, ■ ■ ■ , n. Then, we can consider the chain of 
lattices 

^Aly + ---+Akv ^ ^Aly + ---+An^ C Cy, 
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and we know that £^j^+...+^^^ has dimension equal to k — 1. Hence, we 
have that 

(37) Qu + Q2. + --- + Qk. c /:^,,+...+A. ® K. 

Set Qy := Fo,t, + Qk+iy + • • • + Qnv + Qv We have that 
(38) 

For a sublattice G G C, its saturation will be denoted by s{G). Consider 
the orthogonal decomposition C = s(i3y^i^+...+y^j.^) © C Every point 
p E C may be written as p = q^ +Pkv +P^, where qi, . . . , qind^, are coset 
representatives for £^^^+...+^j.^ in s(£_4^^+...+_4^^), the number of such 
cosets is ind„, pkv e CAi^+-+Ak^ and p^^G £'. 

As (38) is finite, we can intersect Q^ with finitely many rational 
translates of £^^^+...+^^^ IR, and get polytopes 

Q7 ■= (Qmv + Ca,^+...+a,^ ® m) n a, m = 1, ..., M, 

such that (38) is equal to the disjoint union of the sets 

[Qlv + • • • + Qkv + QT) l~l {(imv + Qj + ^Ai^+-+AkJ ) 

where m = 1, . . . , M, andj = 1, • . • , ind^. Now, it is easy to see that 
Q™ := -q,- - qmv + Q™ is a polytope with vertices in £^^^+...+^^^ (g) Q. 
So, as in the previous section, for every p in (38), there exists a unique 
m E {!,..., M} such that p := p — gmu — q?' belongs to 

(39) (Qu + . . . + Qfc. + Qr) n /:^i„+...+A„. 

Consider also the family 

(40) {fi^{Xi,... ,Xn)}l<i<k- 

We may now apply the inductive hypothesis to (39) and construct, as 
in the previous section square matrices M„ q^ indexed by the points in 
(38), whose determinants are non zero multiples of ReSj,(/i„, . . . , /fc„), 
each of them with the same degree in the coefficients of /i^ as the 
resultant. 

As in the primary cell, each of these matrices will have, in the row 
indexed by the point p, the coordinates, in the monomial basis, of 
the expansion of x^'"" fi^{x) for some pair (a, i). In order to define the 
entries of the matrix M, we shall proceed as before: the entry indexed 
by {p,p') wil be the coefficient of x^ in the expansion of x^^"' fi{x). 
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Now that we have the matrix well defined, it remains to decide which 
points lying in secondary cells will be mixed and which not. In order 
to do this, we shall begin with the following 

Definition 4.1. We shall say that the vector v is admissible if there 
exists a unique essential subfamily of (33). 

Recall that Theorem 4.1 in [Stu2] applied to the weight uJ defined in 
(13) gives 



dp 



imtsr(Res^(/o, . . . , /„)) = ± JJ (^^*^j^(/o|fo, • • • , /n|F„)) 

J" 

where JF ranges over all facets of the coherent mixed decomposition 
given hj UJ on Q := Qq + Qi + . . . + Q™, 

J^ = Fo + . . . + F„, 

Fi = conv{A'i}, A'i C Ai, 

ReSjf being the sparse resultant associated with the family of supports 
(^Q, . . . , A'^), and d^ equals the unique integer such that 

(Res^(/o|Fo,--- JuIfJ)'^^ 
has total degree Ya=o MV{Fo, ... , Fi^i, F;+i, ... , F„). 

Remark 4.2. Observe that Res^p- 7^ 1 if and only if JF is associated with 
an admissible vector v. 

The integer d-p will allow us to "choose" the mixed points as follows: 

• if f is not admissible, then all points lying in JF^ fl £ will be non 
mixed. 

• If f is admissible, then we may choose d^^ of the matrices M^ q 
made by subdividing J-'y (1 C in sets of the form (39). We will 
see in Proposition 4.3 that there is at least this number of such 
matrices. 

— If a point indexes one of the d^^ matrices M^ q and is "mixed" 
for that matrix, then it will be mixed for the matrix M. 

— All the other points lying in J-'^Ci C will be non mixed. 

Again, it is easy to see that, in the case n = 1, we get the classical 
Sylvester style matrices ([Syl, Mac]). 

Proposition 4.3. // v is admissible, then the number of sets of the 
form (39) is greater or equal than d-p . 
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Proof of the Proposition: It is easy to see that inituj(Resj[{fo, ■ ■ ■ , fn)) 
is actually the sparse resultant Res^ specialized in the family 

COjbo-^ 1 III ' ' ' 1 In- 

In order to have a nice interpretation of djp , we will use Minimair's for- 
mula given in [Min, Theorem 1] for computing Res^(co,bo2;'''') fir ' ' y fn)- 
First of all, observe that the hypothesis of Minimair's Theorem is sat- 
isfied, due to the fact that the family Aq,--- ,An is essential. This 
implies that the unique essential subfamily of {bo}, Ai, ■ ■ ■ , An is {bo}. 
Using this formula, we have that: 

Res^(co,box''^ /i, ■■■,/«) = Co'Jo n^^s-(/^i^' ■ ■ ■ ' f^J"^ 

V 

where the product ranges over all primitive inward normal vectors of 
the facets of the Newton polytope of ^i + ■ ■ ■ + An, and 

(41) e„ = {apoiv) + apoi-v)) [C^ : Ca^,+-+a^J e^, 

• Fq is the convex hull of {bo} U Aoy, 

• If {^iij,, ■ ■ ■ Afc„} is the unique essential subfamily of {Ai^, ■ ■ ■ , Anv}, 
then e^ is defined as follows: 

1. Let CA,,+-+An. = s{Ca,^^+-+a,^J © C° be the orthogonal 
decomposition, and denote by tt the projection over the second 
factor. 

2. Define e^ := MV {TT{QiXi^{i,,..,i^} , where MV{-) denotes 
the normalized mixed volume with respect to C°. 

Now we will prove the proposition. We can suppose w.l.o.g. that Q = 0, 
because adding a polytope to the Minkowski sum cannot make decrease 
the number of sets of the form (39). Moreover, we may suppose w.l.o.g. 
that the unique essential family is Ai^, ■ ■ ■ , Akv- 

If V is an admissible vector, then e^ and djr^ coincide. Following the 
construction of the Sylvester style matrix for the generalized unmixed 
case (Section 3.1), it turns out that JF^ fl £ may be decomposed as a 
disjoint union of (apoiv) + apoi—v)) [C^ '- '^Ai„+---+A„y] sets parallel to 

Qiv + ■ ■ ■ + Qnv- 

From each of these partitions we may have at least as many subsets 
of the form (39) as integer points are in 

(42) 7r(g,+iJ + --- + 7r(g„J. 

It is well-known that the mixed volume of these polytopes may be 
computed as the number of integer points lying in a subset of a dis- 
placement of (42) (see for instance [Emil, HS]), so e^ will be less or 
equal than the number of integer points in (42). D 
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Remark 4.4. The reader can check that the role played by 6 in the 
previous section is played here by the additional polytope Q : at each 
step of the recursion, we need to impose some conditions on the different 
additional polytopes of lower dimensions in order to guarantee that all 
integer points are in the interior of a cell. It is easy to see that this 
happens for a "generic" polytope. One can set, for instance, 

Q = fixed polytope + 6, 

with 6 generic as in the previous section. 

4.2. Generalized Macaulay style formula. What follows may be 
regarded as the main result of this paper, and an extension of Theorem 
3.8: 

Theorem 4.5. M. is a Sylvester style matrix and det(M) ^ 0. More- 
over, we have the following formula a la Macaulay: 

det (M) = Res^(/o, ...,/„) det (E) , 

where E is the square submatrix of M made by omitting all rows and 
columns indexed by non mixed points. 

As E does not contain coefficients of /o, we get again that det (M) 
has the same degree as Res_4(/o, . . . , /„) in the coefficients of /q. Re- 
placing the rol played by /o with /j, i = 1, ■ ■ ■ ,n, we have a formula 
for computing Res^(/o, . . . , /«) as the gcd of n + 1 determinants. 



Proof. The same argument given in the proof of Theorem 3.8 may be 
applied to this situation in order to see that M is a Sylvester style 
matrix. So, we have that Res^(/o, . . . , /«) divides to det (M) . 

As before, we shall prove that det(M) is not identically zero by show- 
ing that its initial term with respect to ZJ is non zero. Again we have 

degcoeff (/o) (Res^(/o, ••-,/«)) = deg^oeff (/„) (det(M)) . 

Then, 

det(M) _ initjj {det{M)) 

Res^(/o, • • • ,fn) initzj (Res^(/o, . . . , /„)) 

coeffficient of Cq^jj in det(M) 

coefficient of cj^^^ in Res^(/o, ... , /„) 

We will see that the numerator of this fraction is non zero, and that the 
ratio is det(E) by induction on n. The initial case was already covered 
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in the previous section. Suppose n > 1, and let us introduce again a 
parameter of deformation t : 

fi,uj = fi i>l- 

Consider the modified matrix M [ci^a't'^ ) ■ For p E £, let h{p) be the 
largest rational number such that 

{p, h{p)) e Qn = Qo,n + Qi,n + • • • + Qn,n + Qn- 

For every p E S, we shall multiply all the entries in the row indexed 
by p by t^(p)-^i(°-) where -as usual- (i, a) denotes the row content of p. 
We shall denote this matrix with M'(t). It is not hard to see that the 
leading coefficient of det(M'(t)) is initzj{det{M)) . 

The following assertions may be proven mutatis mutandis the results 
given in the previous section: 

Lemma 4.6. It turns out that 

1. < h{p) < 1, for every p E S. 

2. h{p) = 1 if and only if p belongs to the primary cell. 

3. If p and q both belong to the same secondary cell, let us say Tv, 
then 

h{p) = h{q) ^^ {p^v) = {q,v). 

4. If p E J-'v n C has row content (i, a) and v' j^ fiv, /i > 0, then 

(43) {{p - a, h{p)) + Q,,n) n Qniv',i') = 0, V/' G M<o. 

Proposition 4.7. Let < 71 < 72 < . . . < 7Ar = 1 6e the different 
values of h{p) for p G E. Then, the leading coefficient of det(M'(t)) 
(regarded as a polynomial in t) factorizes as follows: 

N 

(44) nd^t(^^-)' 

where Wlj is the square submatrix of M made by choosing all rows and 
columns indexed by those points p such that h{p) = 7^. 
This product may be also factorized as follows: 

(45) det(M^) Yl (det(M;,) . . . det(M^)) , 

:f^ 
donde 

• the product ranges over all secondary cells JF^,, 

• For every v, dy, is the number of sets of the form (39) obtained by 
subdividing T^^ 
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• M^ is a matrix of the type M^g, the latter being defined just after 
(40). 

By the inductive hypothesis, it turns out that det(M^) 7^ 0, Vf, i. 
Using the same argument given in the proof of Theorem 3.8, we also 
get that 

det(M^) = c^^^ det(E^) ^ 0. 

Here, Eat is the square submatrix of Mjv made by all those rows and 
columns indexed by non mixed points. This implies that det(M) 7^ 0. 
Moreover, 

(46) 

mJfedet{M) = c^'^o^dct{Ejv)]^ (dct(Mij . . . dct{M;J^'i )) \[ [dct{ml^) . . .dctimtl'')^, 

where JF^^ ranges over all secondary cells associated to admissible 
vectors f 1, and JF^^ over all cells indexed by non admissible vectors V2- 
Choosing all rows and columns of M indexed by non mixed points, 
a similar version of Proposition 3.12 holds for the matrix E instead of 
M, and it turns out that E will also have a block structure which will 
allow us to compute its determinant as follows: 



(47) 






det(E)=det(E^) J] J]det(E;j J] det (M^ J J] J] det (M^ J . 

From here, the proof of the theorem follows easily, just using Proposi- 
tion 4.3. D 



Corollary 4.8. For every i = 0, 1, . . . , n, the number of i-mixed points 
is exactly MV{Qo, ■■■ , Qi-i, Qi+i, ■■■ , Qn)- 



Corollary 4.9. // the family {.4.j}j=o,...,n is essential, then every co- 
efficient Ci^a appears in Res_4(/o, . . . , /„) with highest power equals to 
■M.V{Qo, . . . , Qi-i, Qi+i, ■ ■ ■ , Qn)- 



Remark 4.10. The same observation given in Remark 3.15 holds also 
here. 
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4.3. Examples. 

Example 4.11. Let us compute the resultant given in Example 2.5. 
In order to do this, we will take bp = (0,0) and Q = {(0, |)}. Here, 

The primary cell is Qi + Q2, and we will have four secondary cells: 



V 


type 


(2,-1) 


2-mixed 


(-1,-2) 


1-mixed 


(-1,-1) 


non mixed 


(-3,-1) 


2-mixed. 



(4J 



In order to subdivide the primary cell, we take uji = (0, 0, 0) and UJ2 = 
(1,1). This lifting produces two cells: a copy of Qi, and the unique 
0-mixed cell of the subdivision. The set S has 23 points, and we will 
associate the unique non mixed cell of (48) with /2. With the aid of 
Maple, we have computed M : 
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Its determinant equals c^Res_4(/o, /i, /2)- The extraneous factor here 
is the first principal minor of size 4x4. 

Example 4.12. Let us compute the resultant of the same system given 
in the previous example, but now we will take Q as the unit square 
[0, 1] X [0, 1] translated by the vector (e, |), with ^ ^ e > 0. The initial 
point bo will be again (0,0). The primary cell will be now equal to 
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Qi + Q2 + Q, and now we will have 6 secondary cells: 



V 


type 


(2,-1) 


2-mixed 


(0,-1) 


non mixed 


(-1,-2) 


1-mixed 


(-1,-1) 


non mixed 


(-3,-1) 


2-mixed 


(-1,0) 


non mixed. 



Observe now that there are secondary cells associated to normal vectors 
of Q which do not define any facet in Qo + Qi + Q2- We subdivide the 
primary cell by taking 



c:'i(0,0) 

C^2(l,l) 

cD(0,0) = 



^^^1(1, 2) 

c:'2(3,o) 

cD(l,0) = 



1, 
0, 



c:'i(2,o) 

Cj{0,1)-- 



= 0, 

4, cD(l,l) 



0. 



If we associate the points lying in the non mixed secondary cells with 
/i, the matrix M obtained has size 36 x 36 and its determinant equals 
cf 6| C2 62 times the resultant. The submatrix E is the following: 
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Example 4.13. We shall see in this example that the inequality estab- 
lished in Proposition 4.3 may be strict. Consider the following essential 
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family: 










Ao = 
Ai -- 
A, - 


= {(0,0);(1,0);(0,1);(1,1)} 
= {(0,0);(1,1)} 
= {((1,0); (0,1)}. 




Denote the 


generic polynomials having those supports 


as follows: 




/o 
/i 


= Oi + a2X + a^y + 04x1/ 
= 61a; + 621/ 






/2 


= ci + C2xy. 





Taking bo = (0, 0), with 1 ^ e > 0, and Q as the triangle with vertices 
(0,0), (1,0), (1, 1) shifted with (e, 1), we get 

• The primary cell is Qi + Q2 + Q- 

• There are five secondary cells, associated with the following in- 
ward vectors: (1,-1), (0,-1), (-1,-1), (-1,0) and (-1,1). 

• (jF(i _!)) nZ^ has two points, while an explicit computation reveals 
that dir =1. 

In this case, {J-'(i-i)) fl Z^ is subdivided in two smaller cells, each of 
them contains exactly one point, and any of the two points may be 
chosen as mixed. If we take Ui, i = 1,2 always equal to 2, uj always 
equal to 0, and apply the algorithm given in the paragraph 4.1, we will 
get the following 13 x 13 matrix: 
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a4 as 03 ai 

We have ordered the rows and columns of M in such a way that the 
monomials lying in (jF(i _i)) fiZ^ index rows and columns 7 and 8. All 
the other non mixed points index the first six rows and columns of M. 
An explicit computation reveals that 



det( 



62C2 6iRes^(/o,/i,/2) 
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and it is straightforward to check that, in this case, h2C^h\ may be 
obtained by taking the principal minor of size 7 x 7, or computing the 
minor indexed by the monomials lying in the first six rows and the 
monomial number eight. 
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